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^ ! Abstract 

A spin-isospin dependent three-dimensional approach has been applied for formulation of the 
three-nucleon bound state and a new expression for Faddeev equation based on three-nucleon free 
basis state has been obtained. Then the three-nucleon wave function has been obtained as a 
function of five independent variables. 
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I. INTRODUCTION 



During the past years, the three-dimensional (3D) approach has been developed for few- 
body bound and scattering problems [lj]-[l3j]. The motivation for developing this approach 
is introducing a direct solution of the integral equations avoiding the very involved angular 
momentum algebra occurring for the permutations, transformations and especially for the 
three-body forces. 

In the case of the three-body bound state the Faddeev equation has been formulated for 
three identical bosons as a function of vector Jacobi momenta, with the specific stress upon 
the magnitudes of the momenta and the angles between them Adding the spin-isospin 
to the 3D formalism was a major additional task which carried out in Ref. jlj]. In this paper 
we have attempted to reformulate the three-nucleon (3N) bound state and have obtained a 
new expression for Faddeev integral equation. To this end we have used 3N free basis state 
for representation of 3N wave function. 

This manuscript is organized as follow. In Sect. [TT] we have derived a new expression 
for Faddeev equation in a realistic 3D scheme as a function of Jacobi momenta vectors and 
the spin-isospin quantum numbers. Then we have chosen suitable coordinate system for 
describing Faddeev component of total 3N wave function as function of five independent 
variables for numerical calculations. Finally in Sect. II III a summary and an outlook have 
been presented. 

II. 3N BOUND STATE IN A 3D MOMENTUM REPRESENTATION 
A. Faddeev equation 

Faddeev eq uation for the 3N bound state with considering pairwise-interactions is de- 



scribed by 



\^) = G tP\^), (1) 



where *) is Faddeev component of the total 3N wave function, M t bing the projection 
of total angular momentum along the quantization axis, P = P12-P23 + -P13-P23 is the sum 
of cyclic and anti-cyclic permutations of three nucleons, t denotes the two-body transition 
operator which is determined by a Lippmann-Schwinger equation and Gq is the free 3N 



propagator which is given by: 

m 4m 

where E is the binding energy of 3N bound state. In order to solve Eq. ([T]) in the momentum 
space we introduce the 3N free basis state in a 3D formalism as [6[: 

|PQ7) = \pqm sl m S2 m S3 m tl m t2 m t3 ) = |qm Sl m tl )|pm S2 m S3 m t2 m t3 ), (3) 

This basis state involves two standard Jacobi momenta p and q which are the relative 
momentum vector in the subsystem and the momentum vector of the spectator with respect 



to the subsystem respectively [14| . (7) = \m Sl m S2 m s:j m tl rn t2 m t:j ) is the spin-isospin parts 
of the basis state where the quantities m Si (m ti ) are the projections of the spin (isospin) of 
each three nucleons along the quantization axis. The introduced basis states are completed 
and normalized as: 

^2 dp dq|pq7)(pq7l = 1, (p'qV|pq7) = <Kp' - p) 5 (q' - q) <W ( 4 ) 
7 J J 

Now we start by inserting the completeness relation twice into Eq. (jTJ) as follow: 

m Am 7" ' 7' 

x(pq 7 |t|p"q'Y)(p /, q , Y|P|p , qV)(p , q , 7 , |^)- (5) 
The matrix elements of the permutation operator P are evaluated as js]: 

(pWl^lp'qV) 

= 8(p" - -q" - q') 6(p + q" + -q') £ m » ml 3 mi 2 

and for the two-body t-matrix we have: 

(Pq7l*|pW> = (P^ 2 ^ S 3 m i2 m t3l t ( e )|p / X / 2 ^3 m tWD^(q _ q // ) 5 m sl m' s ' 1 ^m tl m^, (7) 
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where e = E — -^-q , is the energy carried by a two-body subsystem in a three-nucleon 
system. Substituting Eqs. (JSJ) and (j7]) into Eq. (jSJ) yields: 

(pq 7 | </>»} 

m 4m m' si m tl 

X 



5^ ( P w S2 m S3 m t2 m t3 \ t ( e) | tt m' S2 m' ti m' t2 )(- ir'q[m' ai m' S2 m Sl ro* 1 | ^ Mt ) 

's 2 m' t2 

+ {vm S2 m S3 m t2 m tz \t{e)\ - tt m^m; i m; 3 mi i )(7r / qW Si m Sl m / S3 m; i mi 1 m; 3 |V' Mt ) 



— "3? E / ^ 



£1 £21 

m 4m m' ai m' tl m' s m' t 



x | (p m S2 m S3 m t2 m t3 \t(e)\ir m' si m' s m' tl m' t ) (Tr'q'm' sl m Sl m' s m^m' t m tl \P 2 3\'4> Mt ) 

+ (p m S2 m S3 m t2 m t3 \t(e) P 2 3 1 n rn' si m' s m' tl m' t ) (ix'c{m' si m Sl m' 8 m' tl m tl m' t \ ?p Mt ) j . (8) 

In the last equality we have used the antisymmetry of Faddeev component of the 3N wave 
function as: 

P 23 |^> = -|V> Mt >, (9) 

and also we have considered: 

7r = ^q + q', 7r' = q+-q / . (10) 



The antisymmetrized two-body t- matrix is introduced as [4|: 

a (p'm' S2 m' S3 m[ 2 m[ 3 \t\pm S2 m S3 m t2 m t3 }a = (p' m' S2 m! S3 m! t2 m' t3 \t{l - P 23 )\p m S2 m S3 m t2 m t3 ) , (11) 
where |p m S2 m S3 m t2 mt 3 ) a is the antisymmetrized two-body state which is defined as: 

\pm S2 m S3 m t2 m t3 ) a = -J=(l - P 23 )\p m S2 m S3 m t2 m t3 ) . (12) 
Hence the final expression for Faddeev equation is explicitly written: 

(pqrf\tp Mt ) = — ^ 3^ Yl J d<^ a {vm S2 m S3 m t2 m t3 \t{e)\TT m' s m' s m[m' t ) a 

m 4m m' sl m' t ^m' B m' t 

x {tt'cIw! m si m' s m' t m tl m' t \ip Mt ). (13) 
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As a simplification we rewrite this equation as: 

-1 

E 



1#*(p,q) = - / 5? I d^t a <<^ mt3 (p,n;e)^ «q'), (14) 



m 4m rn'J m' t ' m' a m' t 

where we have used index 7 instead of mgm si m' s rn"mt 1 rn' t for simplicity. This new expression 
is more simple for numerical calculations in comparison with previous expression which has 
been presented in Ref. [5|: 



^a* (P) q) - ~ p 3^ ^ ^ ^T' ' ^ m -3<i ^ 



^1 £2 £<r ' / ■■-. 3 ---s 1 m *3 m t 1 

m 4m 77' a' 



„/ ™/ m / ™/ 



x y rfq'UljI^^^P.-Tri^eiTr'.q'), (15) 

For solving the Eq. (14) one needs the matrix elements of the antisymmetrized two-body 
i-matrix. We connect this quantity to its momentum-helicity representation in appendix A. 
To solve this integral equation numerically, we have to define a suitable coordinate system. 
It is convenient to choose the spin polarization direction parallel to the z axis and express 
the momentum vectors in this coordinate system. With this selection we can write the 
two-body t- matrix and 3N wave function as (see appendices A and B): 

j. m'Jm' g m' t 'm' t / \ 
L a m a2 m S3 mt 2 mt 3 \f' "> fc / 

_ -i[{m S2 +m S3 )^p-{m'J+m' 3 )tp„] f ±m'Jm^m' t 'm' t / rntcn r ir v • f\ (If,) 

~ c L a m S2 m S3 m t2 mt 3 U^> X P> UUB fp^i ^tt, " > i/JOT) \ 1{J J 

^?'(p,q) = e- i [( m ^+^3)^+(^ 1 -^)^ ±Ttff*(p, Xp ,coa(p„,x q ,q), (17) 

< l (ir' ) q') = e^^+W"^ ^fK, a*, cos ^, a/, </)■ (18) 

where x' = q • z, = </V and the labels ± are related to the signs of sin^ p7r , sin<^ pg and 
sin</v<j'- With considering: 

¥V = <f' + Vnq', ¥>'n = V' + ¥n>q>- (19) 

Eq. (14) can be written as: 

± ^ t (p,x p ,cos <p pq ,x q ,q) 

= _ / ^g' / ^ / ^' e «KK)»',,' g-'K^m',)})^,; giKj-AliK}),-})') 

,11111 J-i Jo 

x ^m^jK^s^'^' 008 ^'^'^'^'^ ^^(^^cos^y,^')) (20) 
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where the variables are developed similar to the 3N scattering as 

Xq = q • z, 
Xp = p • z, 



n 



~ = \/^Q 2 + Q' 2 + QQ'Vqq', 

n' = \jq 2 + ^q' 2 + qq'Vqq', 



7T 

3V = 7T • Z = . 



+ q'ypq' 

Vpir = P • TT = — 



7T 

lqyqq' + q' 

Vnqi = 7T • q = , 

7T 

qVqq' + y 

y^q, = 7T • q = . 



Vpq = P • Q = + \/l - A/ 1 - ^ COS <£p<7> 



2/p<3 



COS ip pn 



p ■ q = x p x' + v/l - ^pV 1 - x' 2 cos(y? p - (//), 

Uqq' = q ' q' = X q x' + J 1 - X 2 , V 1 ~ x' 2 COs(y? 9 - V?'), 
p ■ 7T - (p • Z)(7T • Z) ?/ p7r - XpX^ 



n/i- 


(p-z) 2 v 


'l-(*-z)2 ^-4^1-4' 


*' 


■ q' - 


•Z)(q'-Z) _ V^q'-X'^x' 


n/i- 


(tt'-z)S 


/l-(q-z) 2 v/l-x^Vl-x' 2 


7T ■ 


q' - (tt • 


z)(q-z) y^-x^x' 



COSOJw = — — — = = — . 21 

y/1 ~ (TT • Z) VI - (q " Z) 2 V / 1 3 4V / T^ 

It is clear that the Faddeev component of the wave function -0 is explicitly calculated as 
function of five independent variables. In appendices D and E we discuss about the iff- and 
x'-integration and also determination of the signs of sine functions without any ambiguity. 

In this stage we discuss about the total number of coupled integral equations. The total 
number of coupled Faddeev equations for the 3N bound state in a realistic 3D formalism 
according to the spin-isospin states is given by: 

3 

N = 2 (N t x N s ) = 2 (N t x ^ N^), (22) 

i=l 

where N s and N t are the total number of spin and isospin states respectively and N msi is the 
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number of spin states for each nucleon. It is clear that N ms . = 2 and N t = 3 for our problem. 
The factor 2 is related to signs of sine functions of azimuthal angles which is explained in 
appendix D. Consequently the total number of coupled Faddeev equations for either 3 H and 
3 He is iV = 48. 



B. Total wave function 



The total 3N wave function |\l> Mt } is given by [14j : 



\y M *) = (l + P)|V> Mi ). (23) 

Now we derive an expression for the matrix elements of the total 3N wave function by 
inserting the 3N free basis state as fallow: 

<pq 7 |v|> A/ '> = (pq 7 |V Mt > + (pq 1 \P 12 P23^ Mt ) + <pq7l^i3« Mt >- (24) 



By applying the permutation operator P12A3 an d -P13-P23 to the 3N free basis state, Eq. (J2| 
can be written as jf]]: 

(pq 7 |* Mt } = ( P q7l^ Mt ) + (p 2 q 2 72|^ A/t ) + <p 3 <hTs\1> Ut ), (25) 

with: 

1 3 1 
p 2 = --p--q, q 2 = p --q, 72 = m S2 m S3 m Sl m t2 m h m tl , 

1 3 1 
Ps = -2P+4q' q3 = -p-2^' l3 = m S3 m Sl m S2 m t3 m tl m t2 . (26) 

As a simplification Eq. ( |25|) is rewritten as: 

*f (p, q) = <'(p, q) + ifc?(p2, q 2 ) + <<(P3, (27) 
Now we rewrite this equation in the selected coordinate system as: 

*f (p,q) = e -*[(^ 2 +^ 3 )^+(-n-^)^l i^^^COS^X,,?) 

+ e - i[{ms3+ m n ) 9p2 +(m S2 -M t ) Vq2 ) ±^M 4(p2; Xpa , COS X ft , <fc) 

+ e - i [( Wsi+ m S2 )<, J ,3+(m S3 -M t )^3] ±^(p 3 , Xps , COS ^ , X g3 , q 3 ) . (28) 

By considering: 

<Ppa =( Pq + <Ppaq, Vq 2 = ¥q + ¥W 

^P3 =( Pq + Vpsqi fq 3 =¥q + Vqsq, (29) 



Eq. (I28p can be written as: 

± ^* ft (p,x p , cos ip pq ,x q ,q) 
= ± ^(p,x p ,cosi Ppqi x q ,q)+ e^+^K, 

X^e-^+^^+i^-M,)^] ±^M t(p2j Xpa , cos v9 P2(?2 , x, 2 , g 2 ) 

+e _ i[(msi+ms2) ^ 3?+{ma3 _M t) ^ 39 ] ±^t( p3j cos^ m3) x 93 , g 3 )|. (30) 

where: 



cos v? P2g2 



cos ^p 3 g 3 



cos (p P2q 



COS V9 p3 g 



COS (^g 2 g 



1 3 , 1 / „ 9 



= I - -p - -q| = -a/p 2 + -q 2 + 3pqy Pq , 



1 3 , 1 / „ 9 



Pa = I - 2P + = 2 V p2 + ~ 3 ^p«' 



?2 = |p - ^q| = Jp 2 + p 2 - pqy pq , 



qs = I -p-^q| = Jp 2 + \q 2 + pqy Pq , 



X P2 = P2 ' Z 

X P3 = P i ' Z 

= q2 ■ z 

x <?3 = q3 ■ z 



1 3 

2~P x p ^q^q 
P2 

~2~pXp -\- ~^qx q 

P3 

pxp -^q^q 
q2 



93 

P2 ■ q2 - (P2 ■ z )(q2 ■ z ) 



P2?2 ■ L V2 X q2 





(P 2 • 


z )Vi-(q 2 - z ) 2 


Ps 


q 3 - 


- (p 3 ■ z )(qa ■ z ) 


>A- 


(P3 • 


z)Vi-(q 3 - z ) 2 


F>2 


•q- 


(p 2 -z)(q-z) 




(p2 • 


z )Vi-(q- z ) 2 


Ps 


cl- 


(p a ■ z)(q- z) 




ips • 


z )Vi-(q- z ) 2 


q 2 


•q- 


(q 2 -z)(q-z) 


>A- 


(q 2 ■ 


z)Vi-(q- z ) 2 


q 3 


■ q- 


(q 3 -z)(q-z) 



y?\ — Xp 2 ^l — X 



\v 2 -\g 2 -\yqy vq _ 

p 2q2 Xp3 X q 3 



1 — 2^3 a/1 — X 2 



p 2 ^p 2 ^ q 



P3 



g2 ^q 2 - L q 

\A ~ X q2\/^ ~~ x q 

-pypi-h _ 

-i;; m i -■!,! " .< v m g 3 <?3 <? /oi\ 

COS ^39 = /-, ,„ ^ /-, ,„ ^ = 71 2 /l= f ' ( 31 ) 
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The labels ± are related to the signs of sin (p pq , sin (f P2q2 and sin (f P3q3 which are determined 
in appendix D. 



III. SUMMARY AND OUTLOOK 



We extend the recently developed formalism for a new treatment of the Nd scattering in 
three dimensions for the 3N bound state [13]. We propose a new representation of the 3D 
Faddeev equation for the 3N bound state including the spin and isospin degrees of freedom 
in the momentum space. This formalism is based on 3N free basis state. This work provides 
the necessary formalism for the calculation of the 3N bound state observables which is under 
preparation. 



Appendix A. Anti-symmetrized NN t-matrix and its helicity representation 

In our formulation, we need the matrix elements of the anti-symmetrized NN t-matrix. 
We connect these matrix elements to the coresponding ones in the momentum-helicity rep- 
resentation. The antisymmetrized momentum-helicity basis state which is parity eigenstate 
is given by [4j: 

\p;pS 23 \;t 23 }™ = -^(l-P 23 )|p;pS 2 3AMt 2 3> 

= -^(l-^H^+^lpj^aA^W, (A.l) 

Here S23 is the total spin, A is the spin projection along relative momentum of two nucleons, 
£23 is the total isospin and ^23) = \t23T) is the total isospin state of the two nucleons. r 
is the isospin projection along its quantization axis which reveals the total electric charge 
of system. For simplicity r is suppressed since electric charge is conserved. In Eq. (lA.lj) 



P23 is the permutation operator which exchanges the two nucleons labels in all spaces i.e. 
momentum, spin and isospin spaces, and |p; pS^A)^ is parity eigenstate which is given by: 

IP! f>S 2 3\)n = -^(1 + VnPn)\p; P&3A), ( A - 2 ) 
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where is the parity operator, rj n = ±1 are the parity eigenvalues and IplpS^A) is 
momentum- helicity state. The Anti-symmetrized two-body t-matrix is given by |6j: 



™s 2 ™ S3 m t2 m t _ x (AoVp _ A , ) 

t-o m S2 m S3 mt 2 m t;i {P, P j tj — , U( mtri + mto ),fm'_ -4-m' e 



x ^ (1 - r?7r (-) 523+ * 23 ) 

S , 23t23"" 

11 , „,1 1 



x C(o 0^23; m 2 m t3 ) C(--t 2 3;m' t2 m t:i ) 



xC(--5 23 ; m S2 m S3 A ) C(- -S 23 ; m^m'^o) 

x E < 2 a(^) ^ 23 (p, P'; e), (A.3) 

AA' 

where based on momentum-helicity basis states the two-body t-matrix is defined as: 
^v 23 * 23 (p,p';e) = m {p ] pS 23 X-,t 23 \t(e)\p'-,p'S 23 X'-,t 23 r a , 

(AA) 

These two-body t-matrix elements are connected to the solutions of Lippmann-Schwinger 
equation as follow: 

spS 23 iNip , rfSzsf \ JS23 ( r \ 

^S 23 t 23/ I \ Z^AT=-5 2 3 a N\K X PJ U N\'\ X P'J ,nS 23 t 23 / I \ /a r \ 

h\' (p> p ; e ) = - — ^7 — v aa' (p>p, Vpp'i <0. ( A - 5 ) 

where: 



ypp/ = XpXp/ + ^1 - x 2 p ^Jl - x 2 pl cos v^pp'. (A. 6) 

It should be mentioned that the fully off-shell NN t-matrix tlff 3t2S (p,p' ,y vv >\ e), obeys a set 
of coupled Lippmann-Schwinger equations which are solved numerically in Ref. |4j]. Finally 
eq. (A.3) can be written as: 



m 's 2 m 's 3 m 't 2 m 't 3 , / n, _ -i[(m S2 +m S3 )^-(m' +m' 
t-a m S2 m S3 m l2 mt 3 l,P, P , tj — e * 23^ 

Xfc a m B2 m B3 mt 2 mt 3 {P, X p , COS X pp ' , X p i , , €) , ^A. / J 



where the labels ± are related to the sign of sin ip pp t which is determined as: 

sm(p pp > = ±\/l — cos 2 ipppi. (A. 8) 

we consider positive sign for ip pp t G [0,7r] and negative sign for ip pp i e [7r, 27r]. 
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Appendix B. Azimuthal dependency of the 3N wave function 



We introduce the 3N momentum-helicity basis state as: 



|p;p5 , 23 A,q;q5iA) = |p; pS23A)|q; qSiA), 



where: 



S 23 • pIp^A) = A|p5 2 3A), 



S 1 -q|q5 1 A)=A|q 1 S 1 A). 



Thus Faddeev component of the 3N wave function can be written as: 

<'(p,q) = E (pq7|p;P^3A,q;qS 1 A)( P ;p5 2 3A ! q;q5 1 A|^ M ') ) 

S23AS1A 



(B.l) 



(B.2) 



(B.3) 



with considering: 



\f>S 23 \) = R s (p)\zS 23 \)=e- iS ^e- iS ™ 9 *>\zS 23 \), 
|q^A) = R s (q)\zS 1 A)=e- iS '^e-^\xS 1 A), 



(B.4) 
(B.5) 



We have written: 



(pq7|p;pS2 3 A,q;qSiA) = (pq 7 | J R 5 (p) J R 5 (q)|p; zS 23 A, q; z5 x A) 

= {pqi\e- iS ^ p e~ iS ^ e~ iS ^" e^'lp; z,S 23 A, q; z^A) 

= e -im 31 <p q e -i(m S2 +m S3 )^^ pq7 | e -i5| 3 p g-i^jp. %S 23 X, q; Z&A). 

(B.6) 



Also with considering: 



|p; P S 23 A> = Rj p (p)\pz;zS 23 \) = e -^ +s ^ e^+^pz^X), 
|q;q5iA) = Rj „(q)|gz; z^A) = e -^ q +s^ q e -i{L%+s\)9 q ^. 



(B.7) 
(B.8) 



We have written: 



(p;^ 2S A,q;q5 1 A|V^*> 
{pz-iS^qz-zS^Rf^mt^)^) 

(pz; zS 23 A, gz; z^Ale^+W e ^+Sf 3 )^ e ^+sj^ e ^+SfV*|^ t) 

(pz; zS 23 A, gz; z^Ale^+W e ^+sf 3 )^ e <(L 5+ sj 3 ) v , e i(L»+s»)fl, e i(L*+SfV^M t) 

e !M '^(pz; zS 23 A, gz; z^Ale^+^P e i(^+Sf 3 )^ e i(LS+fi?)9,|^Af t ^ (B.9) 
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Consequently Eq. (B.3) can be rewritten as: 

x (p<r/\e^ e *e- iS Z°<\p-,zS 73 \,q;zS 1 A) 

5 23 ASiA 

x (pz; zS 23 \, qz; zSiA|e i(i £ +s ^ e i{L p +s ^^ e i(L v q +sf)e q ^M t ^ (b.10) 
Finally this equation can be written as: 

Vf*(p,q) = e - i[(m ^+ m =3)^+(^i-^W ±^(p,a: p ,cos^ g ,x 9 , g ) (B.ll) 

Appendix C. Parity and time reversal invariance of the total 3N wave function 

In this section we discuss about properties of the total wave function under the parity 
and time reversal invariance. Parity invariance would mean: 

(pqj\V Mt ) = (pq7|P- 1 P.|^ Mt ) = (-p,-q7| J P.|^ Mt ) = (-P,-q7l^ Mt ) 

= (-p, -q 7 |^ Mt ) + (-ft, -q 2 7 2 |^ Mt ) + ("Ps, -q 3 7 3 |^ M >, (CI) 

where we have used P n \^ AIt ) = \ty Mt ) for the 3N total wave function. Eq. (C.I) leads to: 

(p,q 7 |^> = (- p ,-q 7 |^), 
(p 2 ,q 2 7 2 |^ Mt ) = (-p 2 ,-q 2 72|^ Mt ), 

(P3,q 3 73l^) = (-P 3 ,-q 3 73^ Mt )- (C2) 

So we have: 

^(p, x p , cos Vpq , x q , q) = e^ M '~ Mt ^ ± < /t (p, -x p , costp pg , -x q , q), (C.3) 

(p, x p , cos^ p? , x q , q) = e -^- M *)- ±^(p, -x p , cos^ P9 , -x q , q), (C.4) 

where M s = m Sl + m S2 + m Sz . Time reversal invariance of the total wave function can be 



written as 



15|: 



(pq 7 |^ A/t > = (p q7 |T- 1 T|^ Mt ) = 2 2Ms (-p, -q, - 7 |T|^ Mi > = i 2 (^+^*)(-p, -q, - 7 |^- M <). 

(C.5) 
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Considering parity and time reversal invariance lead to: 
(pq 7 |^ M *> = z 2 ^ +M ^(pq,- 7 |^- Mt ) 



(pq, - 7 |r Mt > + (P 2 q 2 , ~l2\r Mt ) + (P 3 q 3 , -73|^- Mt ) • (C.6) 



So we have: 



± ip^ t (p,x p , cos (p pg ,x q ,q) = i 2 ( M *+ M t)* ± tp_^ t ( Pl x In cosip pg ,Xg,q), (C.7) 
(p,x p ,cosip pq ,x q ,q) = iW.+Mt)«±y z Mt (p :Xp ,co S <p pq ,x q ,q). (C.8) 



Appendix D. The (^'-integration 

According to Eq. (120]) the ip'- integration for fixed p, q, x p , x q , cos ip pq and q' can be written 

as: 

/■2tt 

I(ip p ,<p q ) = / dip' e™ 1 ^*-^ e +im2 ^' e- im3 ^'"' 
Jo 

x A ± [cos(ip q - tp'),cos(Lp p - (p'),cos(p pq ]B ± {cos(Lp q - ip')}, (D.l) 



where the and B ± are known functions determined by and 1 * 1 ^ respectively. As we 
know the exponential functions e +ini2LPn i' and e~ m3>Pw 'i' are functions of cos ip nq r and cos ip^ q > 
by considering their sine functions as: 



sin ip nq > = ±J 1 - cos ipl ql , sin yv</' = ±y 1 - cos v . (D . 2) 

Also the cosine functions cosyVg/ and cos p n > q i are function of ip q — ip' . Substituting ip" = 
ip' — ip q leads to: 

1-2-K 

I(ip p , <p q ) = / dip" e~ imiip " e +im2ip ^' e - im ^'*' A* [cos p", cos(^ OT - y>"), cos (p^B^coa ip"] 
Jo 

= I* (cos ippq), (D.3) 

where: 

cos((/9 pg — ip") = cos </? OT cos ip" + sin sin ip" , (D.4) 

and the labels of / ± (cos^ pg ) are depends on the sign of sm.(p pq . It is clear that the angles 
tp nq i and ip n > q > are belong to the interval [— tt, 0] when ip" vary in the interval [0, n] and they 
are belong to the interval [0, 7r] when y/' vary in the interval [k, 2tt]. Furthermore since the 

13 



labels of B^ are depend on the sign of sin </Vy, thus for ip" G [0, tt] and <//' G [vr, 27r] we can 
choose negative and positive labels respectively. Consequently the integral (cos ip pq ) can 
be decomposed as: 

r^cos^) 

/*7T 

= / ^"e" imi ^" e" 1 ™ 21 ^' 1 e +im3| ^'«'U ± [cos(^",cos(^ - (p"), cos ^ M ]S-[cos p"] 
Jo 

/•27T 

+ / efy>"e- <miV " e +im2 *W e- im3 ^' 9 'A ± [cosv? / ',cos(^ - y?"), cos </? M ]5 + [cos y/']. 

(D.5) 

Now we discuss about the labels of A^. As we know the labels of A* 1 are related to the sign 
of sinc^ p7r . We can write p pn = ip pq — (p nq , and then we have: 

sin <p p7T = sin (p pq cos p nq - cos ip pq sin ip nq , (D.6) 

where: 

7T • q- (tt • z)(q- z) y vq - x n x q 



COS (fiirq = 



- (D.7) 

7T 

It is clear that the angle p nq is belong to the interval [0, tt] when ip" vary in the interval [0, tt] 
and is belong to the interval [tt, 2tt] when ip" vary in the interval [tt, 2tt]. Thus depending 
on various intervals of variables p pq and ip", we can choose the positive or negative sign 
for sin<p pq and sinyVq, and then we can calculate simp pn from Eq (B.6). Consequently for 
smtp pn G [0,1] and sinp pw G [—1,0] we can consider positive and negative signs of A^ 
respectively Substituting ip'" = 2tt — ip" , in the second integral of Eq. (B.5) yields: 

pit 

/ dtp'" e +tmiip "' e +im2lp ^' e- ims ^'«'A ± [costp"',cos(ip pq + ip'"), cos ip pq ] B + [cos<p'"], (D.8) 
Jo 

Therefore Eq. (B.5) can be rewritten: 

J^cos^pg) = / dip" e - imiv " e- im2 ^'\ e +im3 ^'o'\A ± [cos ip" ,cos{ip pq - <p"), cos ip pq ]B~ [cos ip"] 
Jo 

+ [ d<p"e +imi<p " e +im2lp ^' e- im ^'«' A^cosip" ,cos{ip pq + iff'), cos ip pq ]B + [cos <//']. 

(D.9) 
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According to Eq. (30) the matrix elements of the total wave function ± \E' can be obtained 
from the matrix elements of Faddeev component of the total wave function ± -0 as follow: 

± ^ 7 [cosv9 p(? ] = ± ^ 7 [cosv? p(? ] + e +imLppq 

x^e" im '^e +im "^ q± ^ 2 [cos^ P2q2 } + e~ m '^ e +in "^" ± ^ 7S [cos <^ 3?3 ] | , 

(D.10) 

It is clear that for (p pq G [0, n], we have: 

ifi P2 q G [7T,27r], tp q2q G [0,7r], 

(p P3q G [7T,27r], ^ G [n,2n], (D-H) 
and for </? P q G [7r, 27r], we have: 

ip P2P G [0,tt], <^ 2g G [tv,2tv}, 

(f P3P G [0,tt], ^ G [0,tt]. (D.12) 

Thus Eq. (D.10) for ip pq G [0,7r] can be written: 

+ ^ 7 [cosy2 P9 ] = + ^ 7 [cosv9 P9 ] + e +imLppq 

x |e +im '^'e w '^ 9 ±^ 72 [ C osv9 P2(?2 ]+ e +m '^ e"™"^ ± ^ 7S [cos </> P3?3 ] j , 

(D.13) 

and for (p pq G [7r, 27r] can be written: 

~^ 7 [cosy? P9 ] = ~^ 7 [cosy? p(? ] + e ~ imippq 

x^e- im '^e- im "^ q± ^ 2 [cos^ P2q2 }+ e~ m '^ e +m "^ q ±V> 73 [cos tp P3q3 ] , j 

(D.14) 

where <pi = 2n — (pi. Now we discuss about the labels of ± ^ 72 and ± V ; 73- As we know the 
labels of ± ^ 72 and ± -0 73 are related to the signs of sin y? P2(?2 and sin tp P3q3 respectively. We 
can write the angles <p P2q2 and (p P3q3 as: 

fp2q2 = fp2q ~~ Pq2qi ^393 = l Pp:iq ~ l Pq:iq- (D.15) 

Consequently we have: 

shiy? P2(?2 = sin^ P2? cos(^ ?2? - cos^ P2 qSin^ ?2? , 

sin ip P3q3 = sin ip P3q cos ip qaq - cos ip Pzq sin ip qsq . (D.16) 
15 



Depending on various intervals of variables tp P2q , f q2q , <£ P3 q and <p q3q we can choose the 
positive or negative signs for their sine functions. Consequently when the calculated sin tp P2q2 
and smcp P3q3 from Eq. (D.16) are belong to interval [0, 1] we can use positive sign of ± -0 72 
and ±- 73 and when they are belong to interval [—1,0] we can use negative sign of them. 

Appendix A: Appendix E. The ^'-integration 

According to eq. ([20]) the x'-integration carried out as: 

± ip(p, x p , cos ip pq , x q , q) = J dx'C(x') ^(tt' ,x^, cos cp^q',^ ,q'), (E.l) 



where the C is known function determined by and exponential functions. This equation 



a iiLin mil! unrjiiiujni in ' 

can be rewritten as: 



~-?p(p,x p , cos ip pq ,x q ,q) = J dx'C(x') ± if; (n'.x^, cos ip^q^x'.q 1 ) 

+ / dx' C(x') ± i(j (it' , x n ' , cos (fij./g/, x' ,q') 
'0 

dx' C{— x') ± ijj(n'(—x'),x n i(—x'), cos (p^ q i(—x'), —x', q ) 

+ i dx 1 C(x') ip(7r', av, cos yv 9 '> x> ■, <l')i (E-2) 
Jo 

Finally by considering parity invariance which is described in appendix C, Eq. (E.2) can be 
written: 



ip(p,x p , cos Lp pq ,x q ,q) 
dx' < (— ) s+ *C(— 2/) ± ifj(ir'(—x'),—x. K '(—x'),cosip 7T i q i(—x'),x',q) 







+ C(x') ± -0(vr',x 7r /,cos(/? 7r ' g /,a; / ,g / )|. (E.3) 
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